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1. Introduction 

The generation of optical frequency combs using mi¬ 
croresonator devices is currently a hot topic that is at¬ 
tracting a significant amount of research. A mode-locked 
optical frequency comb can combine an ultra broadband 
spectrum with coherent and finely resolved comb lines at 
an equidistant frequency spacing. Microresonator based 
optical frequency combs have a plethora of potential ap¬ 
plications in, e.g., the areas of precision metrology, spec¬ 
troscopy, optical clocks, wavelength-division multiplex¬ 
ing sources, and arbitrary waveform generation 0- By 
enclosing a nonlinear medium within a resonator struc¬ 
ture, a wealth of dynamical behaviours appear. Mi¬ 
croresonators with a Kerr nonlinearity that are pumped 
by a continuous wave (CW) laser are known to support 
bistable CW solutions, as well as both periodic tempo¬ 
ral (Turing) patterns and localised cavity soliton solu¬ 
tions. These solutions may be either stable or unstable, 
depending on the operating regime determined by the 
power and detuning of the pump laser. Additionally, 
there exist operating regimes where the intracavity field 
is chaotic. 

Previous work on the generation of optical frequency 
combs in microresonator devices has primarily focused 
on the regime where the intracavity power is relatively 
low. In this case the mean-field approximation is valid, 
and the dynamics can be accurately modeled using ei¬ 
ther the Lugiato-Lefever equation (LLE) 0 or using the 
theory of coupled mode equations 0, 0 . For these ap¬ 
proaches to be applicable, it is assumed that the over¬ 
all cavity detuning, including both linear and nonlinear 
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contributions, remains small during propagation. Par¬ 
ticularly, it is assumed that the characteristic nonlinear 
length scale is much longer than the path length of the 
cavity. In the present work, we aim to go beyond the 
LLE and the mean-field approximation, and consider 
the resonator dynamics when the nonlinear phase-shift 
is relatively large (of the order of unity). 

As a motivation for why it might be prudent to con¬ 
sider microresonator dynamics beyond the LLE, we note 
that mean-field models have been unable to account for 
the complete range of dynamical behaviours that have 
been experimentally observed in resonator devices. Re¬ 
cent experimental work at NIST has e.g. shown exam¬ 
ples of frequency combs that are not well described by 
the LLE Q. Mean-field models are also unable to ac¬ 
count for the period doubling behaviour, which has been 
observed to occur in dissipative fiber-ring resonators 0- 
@]. This period doubling behaviour is distinct from the 
periodic breathing of cavity soliton type solutions of the 
LLE, which occurs on a different time-scale that is not 
associated with the roundtrip time. Period doubling in 
fiber-ring cavities has been experimentally observed in 
the periodic switching of the amplitude of output pulses 
when a fiber-ring resonator is pumped by a synchronised 
pulse train with large input amplitude [§]. The period 
doubling sequence is also well known to be an essential 
component in the route to chaos 0E3. 

While the lengths of microresonators are generally 
about four or so orders of magnitude shorter than fiber¬ 
ring resonators, they can display nonlinearities that are 
several orders of magnitude greater than the x ^ nonlin¬ 
earity of silica fiber. This is especially true for microring 
resonators made of semiconductor materials. The use of 
highly nonlinear microring resonators with fairly long 
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path lengths, in combination with intense input pump 
powers, could consequently enable the observation of mi¬ 
croresonator dynamics similar to what has previously 
been observed in fiber-ring resonators. 

In the next section [2] we introduce the Ikeda map, 
which is the basic model system used to describe fre¬ 
quency comb generation beyond the LLE. We use the 
Ikeda map to analyse the CW solution in section [3] and 
show that additional stationary states with higher intra¬ 
cavity power appear when compared to the LLE. The 
stability of these stationary states is investigated in sec¬ 
tion QJ showing the presence of new parametric instabili¬ 
ties other than the conventional modulational instability 
for high pump powers. In Section [5] we use the para¬ 
metric instability to show examples of period doubling 
behaviour. Thereafter, in section [Gl we demonstrate the 
intriguing possibility of forming super energetic cavity 
soliton associated with the excited states. Finally, in 
the last section 0 we summarise the article and present 
some overall conclusions. 

2. The Ikeda map 

In this section we introduce the Ikeda map 0 , which 
we use to model the dynamics of optical frequency combs 
without resorting to the mean-field approximation. The 
map is constructed by combining an ordinary nonlinear 
Schrodinger equation, for describing the evolution of the 
envelope of the electric field within a waveguide, together 
with boundary conditions that relate the fields between 
successive roundtrips and the input pump field nm, 
viz. 

E m+1 {t, 0) = VeE in + VT L), (1) 

dE m (t,z) a ir , mu , 

dz 2 1 ’ ] 2 dt . 2 

+ i 1 \E m {t,z)\ 2 E m {t,z). (2) 

Here the independent variables are the evolution vari¬ 
able z, which is the longitudinal coordinate measured 
along the waveguide and t which is the (ordinary) time. 
The first equation is the boundary condition that deter¬ 
mines the intracavity field E m+1 (t, z = 0) at the input of 
roundtrip to + 1 in terms of the field from the end of the 
previous roundtrip E m (t, z = L) and the pump field Ei n . 
The path length of the resonator is assumed to be equal 
to L. Additionally, 6 is the coupling transmission coeffi¬ 
cient and </>o = 2irl — Sq is the linear phase-shift, with 5q 
the frequency detuning from the cavity resonance clos¬ 
est to the pump frequency, assumed to correspond to 
the longitudinal mode number 1 = 0. Eq. m is writ¬ 
ten in the reference frame moving at the group velocity, 
with second-order group velocity dispersion coefficient 
and nonlinear coefficient 7 = n 2 uJo / {cA e f f) (with 712 
being the material Kerr coefficient, ui 0 the angular pump 
frequency, c the speed of light in vacuum and A e f / the 
effective mode area), as well a linear absorption coef¬ 
ficient cnj. In this work we will for simplicity neglect 
higher-order dispersion and any frequency dependence 
of the nonlinear coefficient. However, since Eq. m is 


simply the evolution equation of a straight waveguide 
it can trivially be extended to include any higher-order 
effects. 

In the usual treatment of microresonator frequency 
combs, cf. Si, one proceeds by averaging the map 
Eqs.dlK]) over one roundtrip in order to obtain the LLE 


. dE(t,r) ,/3 2 Ld 2 E(t 1 r) 

R dr +l 2 dt 2 

- (a + iS 0 )E(t , r) + V6E in , 


i"/L\E(t,T)\ 2 E(t, t) = 

( 3 ) 


where the evolution variable z has been replaced by the 
slow time r. The coefficient ta is the roundtrip time and 
a = (a.iL + 9 )/2 is the total cavity loss. 

Since the LLE can be derived from the Ikeda map, 
the latter is consequently more general. The Ikeda map 
is easy to simulate using conventional numerical meth¬ 
ods for the nonlinear Schrodinger equation such as the 
split-step Fourier method }l.3| j. However, numerical sim¬ 
ulations of the Ikeda map must use a step length that 
is shorter than the cavity length, which makes it more 
computationally demanding than the LLE. 

The derivation of the mean-field approximation is 
valid provided that the total cavity detuning is much less 
than unity, i.e. 7 LI-E ) 2 -C 1 and S 0 -C 1, see Ref. 0 . 
The latter condition signifies that the pump detuning 
from resonance is small, while the former states that 
the nonlinear phase-shift is also a small quantity, which 
can be equivalently interpreted as a requirement for 
the nonlinear length scale L n i = 1 /( 7 |K| 2 ) to be much 
longer than the cavity circumference. The average dis¬ 
persion over each roundtrip should additionally remain 
small. Note, however, that this does not necessarily im¬ 
ply that the field must be changing slowly inside the 
cavity, cf. [lj. Below we will investigate some of the 
new dynamical regimes that appear as the nonlinear 
length scale becomes comparable to the cavity length, 
i.e. Lni ~ T. 

3. Multi-valued stationary solutions 

In this section we consider continuous wave solutions of 
the Ikeda map Eq. m • Particularly, we will consider 
CW solutions where the intracavity field can decay along 
the resonator circumference due to intrinsic absorption 
(cq > 0 ), but where the same amplitude is reproduced 
after each roundtrip with the coherent addition of the 
pump field. Fields that are periodically restored can be 
considered stationary in the averaged sense of the mean- 
field theory. 

The CW solution of the nonlinear Schrodinger equa¬ 
tion © that accounts for the field that has propa¬ 
gated around the resonator is given by E m (t,L) = 
E m (0)exp(-a i L/2+i 7 L e f f \E m (0)\ 2 ),withL eff = (1- 
e ~ aiL )/ai being the effective nonlinear length due to in¬ 
trinsic absorption (T e // —> L as a, —> 0). Assuming that 
E m+1 ( 0) = E m (0) = E we find that the intracavity field 
should satisfy the following fixed point equation 

E = pe^E + V0E in 


( 4 ) 
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which has a simple geometrical interpretation, see Fig. 
□ and Ref. [T3|. Meanwhile, the intracavity power must 
satisfy the following condition 

1 = JEL (5) 

(1 - p) 2 + 4psin 2 (<^/2) 9\E in \ 2 
0 = 6 0 - r yL e ff\E\ 2 

where the roundtrip loss coefficient p = 
\/l — Qe~ aiL ! 2 k, 1 — a. The right-hand-side of 
Eq. © has a linear dependence on the intracavity power 
\E\ 2 while the left hand side can be recognised as an 
Airy function describing resonances whenever 4> = 2nn 
with n an integer. 

!m(E) 



Fig. 1. Geometric illustration of a fixed point of Eq.©. The 
original vector is reproduced after applying a rotation, a scal¬ 
ing and adding the pump vector. Either one or three sta¬ 
tionary states are possibile when the rotation (phase-shift) is 
assumed small, but additional stationary states can also be 
obtained if the rotation is by more than one lap. 


Log 0|E in | 2 



Fig. 2. Comparison of pump dependence on the intracavity 
power for multi-valued stationary continuous wave solutions 
of the LLE (blue dashed) and the Ikeda map (orange solid). 
Parameters a = 9 = 0.1 and 7 = L = Pi n = So = 1. 

The corresponding bistable steady-state condition for 
the LLE (J3]) is given by 


a 2 + (So - 7L\E\ 2 ) 2 9\E in \ 2 


A comparison of Eqs.©© is shown in Fig. [2j where it 
is seen that the LLE (blue dashed curve) is an excellent 
approximation for small intracavity powers but that it 
quickly becomes bad as the power increases. Indeed, 
from Eq. © we see that the intracavity power will have 
a linear asymptotic dependence on \E\ 2 , while Eq.© 
for the LLE predicts that the intracavity power should 
depend asymptotically on the cubic root of the pump 
power. The intensity inside the resonator will conse¬ 
quently increase much faster than the LLE predicts. 
From Fig.[5]we also see that higher-order multi-stable (or 
excited) states do not necessarily imply that the pump 
power needs to be large. In fact one finds that excited 
solutions may even coexist with the low-power bistable 
response. 

The successive intracavity powers for which the pump 
power has a minimum in Fig. [2j i.e., the least pump 
power needed to reach a state on a given level, can be 
found explicitly from Eq.© as \E\ 2 = (<5 0 + 27 m)/( 7 L), 
where n is a positive integer, while the corresponding 
pump powers are given by 0\Ei n \ 2 = (<5o + 27rn)(l — 
p) 2 /(^L). The pump power required to reach the first 
excited level is thus a fraction 1 + 27 t/ <5o larger than 
that required in order to reach the upper level of the 
bistability curve. Therefore, the observation of excited 
states is facilitated if the detuning is fairly large. In this 
article we refer to states with n > 1 as excited states, 
since the LLE is only able to model the first minimum 
occurring for n = 0 . 

We now consider the physical parameters necessary to 
obtain a nonlinear phase-shift of order unity in microres¬ 
onator devices. The phase condition (t>NL = 7 LP 0 ~ 1 
depends on the intracavity power Pq = \E\ 2 rather than 
the pump power. However, from Eq.© one can esti¬ 
mate that the power will grow linearly with the res¬ 
onator finesse, i.e. Po ~ PinSF / 7 r. This approximation 
corresponds to a lower bound for a critically coupled res¬ 
onator. The phase-shift condition can alternatively be 
written in terms of the resonator Q-factor and the effec¬ 
tive mode area as 6nl ~ 2— ~ 1. The nonlinear 

phase-shift is related to the Kerr frequency mode shift 
as Aitvl = 4>nl^vfsr/( 27t), with a phase-shift of 27 t 
thus able to shift the field to a different resonance. Note 
that the modes can also experience a thermally induced 
mode shift, however, this is not taken into account in 
the current treatment. 

Microresonators featuring large mode shifts have in 
fact already been experimentally demonstrated in the 
literature. For example, in Ref. [© a fused silica mi¬ 
crotoroid, with a diameter of 80 pm and Q « 2.7 x 10 s , 
was used to demonstrate octave spanning comb gener¬ 
ation. This resonator had a Kerr frequency shift that 
was estimated to be as much as 100 GHz at 1 IK of in¬ 
put power, while the FSR of the resonator was 850 GHz. 
The corresponding nonlinear phase-shift would therefore 
be (j)NL ~ 0.74, which suggests that the mean-field ap¬ 
proach may be inaccurate. The comb bandwidth of the 
same resonator was also theoretically analysed in Ref. 
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fl6l |. with experimental measurements showing devia¬ 
tions of one order of magnitude from predictions based 
on the LLE. These large deviations were however at¬ 
tributed to be mainly due to higher-order dispersion. 


4. Linear stability analysis 

Optical frequency comb generation in microresonators 
pumped by a CW laser is initiated by the growth of 
frequency sidebands above a certain threshold power. 
These sidebands are symmetrically spaced and occur 
due to modulational instability (MI) of the pump mode. 
Because of the cavity boundary conditions, MI in a res¬ 
onator cavity can display qualitatively different features, 
compared to the conventional MI in a straight waveg¬ 
uide. This is because the detuning between the reso¬ 
nance frequency of the cavity and the pump frequency, 
introduces an additional degree of freedom for the phase¬ 
matching of different four-wave-mixing processes. As a 
consequence, the steady-state CW solution can become 
modulationally unstable also in the normal dispersion 

regime E3- 

An analysis of the linear stability of the steady-state 
solution in the mean-field limit has been performed by 
many authors, see e.g. Refs. [!, EH, HH . However, the 
LLE does not model the full dynamics of the instabili¬ 
ties that can occur in a resonator cavity described by the 
Ikeda map. Particularly, one finds that additional unsta¬ 
ble frequency bands (resonance tongues) appear at high 
power. These bands can be associated both with insta¬ 
bilities having the fundamental period and with period 
doubling instabilities [20]. 

In order to investigate the linear stability of steady- 
state solutions described by Eq.0, we assume an 
ansatz of the form E m {t,z) = [Eq + u m (t,z ) + 
iv m (t,z )] exp[-a i z/2+i'y(l-e~ aiZ )\E 0 \ 2 /a i +iArg E 0 \, 
where the perturbations u m (t,z ) and v m (t,z ) are real 
functions. Linearising the nonlinear Schrodinger Eq.© 
we find that the Fourier transform of the perturbations 
w m (z ) = [u m (w, z), v m (uj, z)] T should satisfy the follow¬ 
ing system of linear equations dw m /dz = Mw m , with 
the perturbation matrix 


M{z) 


0 

(/3 2 /2 )w 2 + 2 7 |E 0 | 2 e-“^ 


-(/V2V 2 

0 


• ( 7 ) 


The analytical solution of this system is quite involved 
in the general case and the stability is most easily inves¬ 
tigated numerically. However, in the absence of intrinsic 
absorption, i.e. ai —> 0, the matrix M becomes indepen¬ 
dent of z and the solution simplifies considerably, see 
Refs. Himi- In this case one finds that the eigenval¬ 
ues of M give the conventional MI gain of a straight 
waveguide, viz. 


li = -p2l\E 0 \ 2 - (/3|/4)w 2 , 

while the solution of the linear system is given by 

—/?2 w 2 e /iz —/32W 2 e _MZ 
2^ 2 —2^e~ fJ,z 


( 8 ) 


w m (z) = 



Q"m 


bm 


= W(z)c n 


( 9 ) 


In the next step we consider the boundary condition 
Eq.©) that the perturbed solution must satisfy and ap¬ 
ply the Fourier transform. This leads to w m+1 (0) = 
\/l — 0R{^>)w m {L), where <f> = <j>o + 7 L|Eo| 2 and R{<f>) 
is the rotation matrix 


R(<P) 


cos </> — sin (f> 
sin </> cos cf) 


( 10 ) 


We then obtain a system of difference equations for 
the coefficients c m = \a m , b m ] T , viz. c m+1 = 
\/T^0[W- 1 (O)R(cfi')W(L)\c m = Qc m . The stability of 
solutions to this linear system of difference equations 
with constant coefficients depends, similarly to Eq. 0, 
on the magnitude of the eigenvalues. The Ikeda map has 
an instability whenever the modulus of an eigenvalue of 
Q is larger than unity. 

For the case when the intrinsic absorption is neglected 
these eigenvalues are explicitly given by 

q± = V 1 - 0(p± i/p 2 - 1) (11) 

where p = cos(^) cosh(^L) 

[yl-Eol 2 + (/? 2 / 2 )w 2 ] sin(</>) sinh(^L)//x, with p given 
by Eg.© and the phase-shift (j> = </> 0 + 7 L|Eo| 2 , see 
Ref. j20|E 

For the general case the stability of the Ikeda map 
can be determined using a numerical Floquet anal¬ 
ysis E3- This is accomplished by investigating 
the stability of eigenvalues of the fundamental matrix 
[u>™ +1 (0), w™ +1 (0)]. This matrix is obtained by nu¬ 
merically integrating the perturbed system dw m /dz = 
M(z)w m over one roundtrip, starting from two lin¬ 
early independent initial conditions, e.g. w™ 2 {0) = 
[1, 0] T , [0, 1] T , and applying the boundary condition 
Eq.0. 

In Figs. [3J4] we show some some examples of modu¬ 
lational instability gain for both anomalous and normal 
dispersion. The gain has been calculated numerically 
using the Floquet method for parameters corresponding 
to the silica resonator of Ref. [15( mentioned at the end 
of section [3] Note that even though the magnitude of 
the intracavity power is very large, the maximum plot¬ 
ted value may actually correspond to pump powers as 
low as 2 W, owing to the very high finesse T ~ 10 6 of 
the resonator. Fig. [3] shows that the modulational in¬ 
stability analysis for the LLE is valid for low intracavity 
power, but that new unstable sidebands appear as the 
power increases. The high power parametric instability 
tongues may correspond to instabilities occurring for ei¬ 
ther resonant or anti-resonant conditions. The former 
includes the ordinary CW-MI of the LLE and higher or¬ 
der sidebands, and has its origin in the q + eigenvalue, 
while the latter comes from the eigenvalue and gives 
rise to a period doubling instability referred to as P2-MI, 
see Refs. 0, [23}- The P2-MI can be used to convert a 
continuous wave solution to a period doubled pattern, 
i.e. an alternation between two different modulated tem¬ 
poral patterns having a period of one roundtrip each. 
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Fig. 3. Parametric instability tongues of the Ikeda map 
for anomalous dispersion. The red contour shows the pre¬ 
dicted range of modulational instability for the LLE. Below a 
cross section corresponding to the dashed line is shown with 
alternating CW-MI/P2-MI gain bands. Parameters: 7 = 
25 W~ 1 km~ 1 , L = 2n x 40 pm, T = 10 6 , = — 40 ps 2 km~ 1 
and = 0 . 
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Fig. 4. Parametric instability tongues of the Ikeda map for 
the same parameters as in Fig. [3]but with the sign of the dis¬ 
persion reversed to give normal dispersion $2 = 40 ps 2 km~ 1 . 
The cross section shows alternating P2-MI/CW-MI gain 
bands. Note that the LLE does not predict any instability 
for these parameters. 


Since an optical spectrum analyser perforins an averag¬ 
ing over many roundtrips, it may be difficult to exper¬ 
imentally distinguish a stable period doubled frequency 
combs from an ordinary frequency comb. 

While the LLE gives a good approximation for the 
modulational instability gain occurring at low powers, 
its accuracy becomes increasingly bad as the power in¬ 
creases. Indeed, the LLE is not able to model either 
P2-MI instabilities nor higher order CW-MI sidebands. 
Fig.[4]shows that the period doubling instability is in fact 
usually the first instability to occur for normal disper¬ 
sion cavities, except for the case when the detuning has 


been used to compensate for the phase-mismatch. More¬ 
over, we see that high intensity continuous wave solution 
may actually be stable in certain narrow regions. 

The excited continuous wave steady-state solutions 
generally correspond to the parameter regions which are 
stable to plane wave perturbations. However, the stabil¬ 
ity analysis shows that the new states are usually rnodu- 
lationally unstable (CW-MI and/or P2-MI) for anoma¬ 
lous dispersion, similar to the upper steady-state level in 
the simple bistable case of the LLE. But since the mod¬ 
ulational instability only applies to structures that are 
sufficiently broad in time, it suggests the intriguing pos¬ 
sibility of forming cavity solitons associated with these 
excited multi-stable states. Indeed, this turns out to be 
possible under appropriate circumstances and we shall 
see an example of such excited cavity solitons in section 

E3 

5. Period doubled patterns 

In this section we give some examples of period doubled 
patterns. A period doubled CW solution should satisfy 
a system of equations similar to that of a steady-state 
fixed point solution, viz. 

A = pe i * B B + V6E in , B = P e irt,A A + VeE in . (12) 

Here we have defined the half-period fields as A = E m 
and B = E m+1 . The presence of a period doubled fixed 
point implies that the field satisfies a functional relation 
of the form x = F(F(x)). If the fields in each half-period 
are equal, i.e. A = B = E (or x = F(x )), there is no 
period doubling and we recover Eq.dH). From Ea. (fT2l) 
one finds that the two fields are related to each other 
through the following resonance condition, see Ref. [24j , 


1 + pe^ B 1 + pe^A ’ v ' 

where 4>a,b denotes the phase-shift with the nonlinear 
term evaluated for A or B , respectively. 

Period doubled states can occur for relatively low in¬ 
tracavity powers \A\ 2 and \B\ 2 if the cavity detuning is 
close to anti-resonant conditions, i.e. (j> s=s n. The above 
condition m also shows that period doubling cannot 
occur in a linear medium since the denominators cancel 
each other in the absence of a nonlinear phase-shift. 

Two coupled conditions for the intracavity power can 
be derived, viz. 


e\E in \ 2 = \A\ 2 


(1 - p 2 ) 2 + 4p 2 sin 2 l((j) A + <t> B )! 2 ] 
(1 + p) 2 - 4psin 2 (</> B /2) 


(14) 


for A and an analogous condition for B which is obtained 
by interchanging A and B in the above equation. After 
solving these equations, the fields can be recovered using 
the relation A = \/0E in ( 1 + pe i<t>B )/{ 1 - pV^+^s)). 

Period doubled CW solutions are usually only ob¬ 
tained for large input pump powers in microresonator 
systems where the finesse is high. However, other pe¬ 
riod doubled patterns can be found at rather modest 
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power levels as a result of period doubling instabilities. 
An example is shown in Fig. [5] for a microresonator with 
normal cavity dispersion. The period doubled pattern 
consists of two identical but out of phase patterns that 
alternate between each roundtrip. The pattern is found 
to be stable for the considered operating parameters, but 
it becomes unstable due to the growth of additional side¬ 
bands if the pump power is increased further. Note that 
the spectrum is constant for the current example since 
only the spectral phase changes between each roundtrip. 



Fig. 5. Period doubled pattern obtained through P2-MI 
instability for normal dispersion resonator. Two identical 
but out of phase patterns alternate between each roundtrip 
(blue solid line denotes even numbered roundtrips, red dotted 
odd). Numerical simulation of the Ikeda map using parame¬ 
ters a = 6 = 0.01, 7 = 100 W/m, L = 0.1 mm, Pi n = 4.5 W , 
FSR = 100 GHz, /?2 = 500 ps 2 /km and So = 0.25. 


6. Localised solutions, cavity solitons 

An important question that we have not considered so 
far is whether or not the Ikeda map has any additional 
types of localised solutions different from those of the 
LLE. That the answer to this question should be in the 
positive can be appreciated from the following line of 
argument. 

It is well known that cavity solitons exist for almost 
the same range of parameters as the bistable CW solu¬ 
tion in microresonators described by the LLE. The rea¬ 
son for this is that a bright cavity soliton consists of a 
solitary pulse sitting on top of a constant background. 
This background is a stationary and modulationally sta¬ 
ble solution of Eq.(0. Similarly, since cavity solitons 
are concave at their peak, there should also exist an¬ 
other stationary CW solution that corresponds to two 
points on the soliton were the second derivative is zero. 
However, it should be noted that this is only approx¬ 
imately true, due to the phase difference between the 
soliton and the background. The two points are inflex¬ 
ion points, i.e. points on the soliton profile were the slope 


of the derivative goes from increasing to decreasing or 
vice versa. The latter stationary CW level is however 
not required to be stable against MI as long as the soli- 
ton is not too broad. Consequently, it can be understood 
that one should be in a regime were at least two CW so¬ 
lutions coexist for the same pump parameters in order 
to find a stationary cavity soliton solution. Similarly, if 
more than two stationary CW states are present simul¬ 
taneously, then we can expect that each new level may 
also be associated with a different cavity soliton state. 

To find a cavity soliton associated with an excited 
CW state we consider operating parameters allowing for 
three simultaneous steady-state CW levels with positive 
slope, see Fig. [6] For the numerical simulations we con¬ 
sider a microresonator with parameters a = 0 = 0.01, 
7 = 100 W/m, L = 0.1 mm, FSR = 100 GHz, 
fa = —500 ps 2 /km, Pi n = 1 W and <5o = 0.6. These pa¬ 
rameters are somewhat idealised since we neglect higher 
order dispersion, but similar parameters could be ob¬ 
tained for semiconductor microrings featuring large non¬ 
linear Kerr coefficients such as e.g. AlGaAs. In Figs. OH] 
we demonstrate a case were two different stable cavity 
solitons, one conventional LLE soliton and one super 
cavity soliton associated with an excited state, are ob¬ 
tained and may coexist simultaneously inside the cavity. 

The two stationary solitons in Figs. OH] are found in the 


Log S|Ei„| 2 



Fig. 6 . Configuration for stationary CW solutions allowing 
for multi-level cavity solitons. Logarithm of pump power 
plotted vs. intracavity power. The constant line corresponds 
to the external pump power (1 W) used in the examples. 

same microring using identical pump parameters. The 
second cavity soliton that is associated with the excited 
CW level around \E\ 2 ss 68.8 W in Fig. [6] is seen to 
be much more energetic than its LLE counterpart, that 
is associated with the CW power level \E\ 2 ss 6.4 W, 
corresponding to the upper bistable CW solution of the 
LLE. The soliton shown in Fig. [5] has a peak power that 
is one order of magnitude greater than that of the LLE 
soliton. At the same time it is also narrower, leading 
to a mode-locked frequency comb spectrum that is sig¬ 
nificantly broader than that of the LLE soliton. This 
is particularly interesting for applications to broadband 
frequency comb generation, since such super cavity soli- 
tons could easily span a bandwidth of more than an oc¬ 
tave, and possibly even multiple octaves with appropri- 
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Fig. 7. LLE type cavity soliton with background level \E\ 2 ss 
0.029 W and support level |i ?| 2 « 6.365 W. Simulation of 
Ikeda map Eqs. CEE}, with parameters a = 0 = 0 . 01 , 7 = 
100 W/m, L = 0.1 mm, P in = 1 W, FSR = 100 GHz, 
P 2 — —500 ps 2 /km and So = 0.6. 



Fig. 8 . Stationary super cavity soliton with background level 
|i ?| 2 ss 0.029 W and support level |.E | 2 « 68.764 W. using 
the exact same parameters as in Fig. [7] except for different 
initial conditions. 


ate dispersion engineering. Other potential applications 
include multi-state optical buffer memories, cf. [25j . 

The intensity of each soliton sans background can be 
approximated as \E(t)\ 2 = Posech 2 (t/To), with peak in¬ 
tensity Pq = 2(5o + 2nn)/( r yL) and temporal duration 
Tq = (2(<5o + 27m)) for n = 0,1. Note that the 

two solitons were excited numerically by starting with 
initial conditions close to the final soliton profile. This 
corresponds experimentally to the injection of a soliton 
pulse into the cavity. At present it is not clear how super 
cavity solitons could best be excited in an experimental 
setting. This and the relation of super cavity solitons to 
other types of bistable solitons, cf. [261 ]. will be the focus 


of future work. 

7. Conclusions 

In this article we have presented an analysis of frequency 
comb generation that goes beyond the mean-field ap¬ 
proximation used in the Lugiato-Lefever equation and 
coupled mode theory. The analysis has been carried out 
in the context of microresonators but is valid also for dis¬ 
persive fiber-ring resonators. The LLE is valid only in 
the limit of small phase-shifts were the nonlinear length 
scale is much longer than the length of the cavity. Us¬ 
ing the more general Ikeda map as a model, we have 
shown that new multi-valued stationary states appear 
at high intracavity powers for the CW solution. These 
states may be present simultaneously as the low power 
bistable response and do not necessarily require a large 
pump intensity. The modulational stability of the sta¬ 
tionary CW states has been analysed, while taking into 
account the boundary conditions imposed by the cav¬ 
ity. The stability analysis shows the presence of several 
parametric instability bands, including period doubling 
instabilities, at high intracavity power both for normal 
and anomalous dispersion. An example of how such pe¬ 
riod doubling instabilities can lead to a stable pattern 
that alternates between each roundtrip has been demon¬ 
strated. Finally, we have predicted a novel type of super 
cavity soliton solution associated with an excited CW 
steady-state. These solitons are narrower and more en¬ 
ergetic than those of the LLE, allowing for efficient gen¬ 
eration of mode-locked ultra wideband frequency combs. 
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